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^ . Abstract 

, We investigate toy dynamical models of energy-level repulsion in quantum eigenvalue se- 

quences. We focus on parametric (with respect to a running coupling or "complexity" parameter) 
stochastic processes that are capable of relaxing towards a stationary regime (e. g. equilibrium, 
. invariant asymptotic measure) . In view of ergodic property, that makes them appropriate for the 

^ ' study of short-range fluctuations in any disordered, randomly-looking spectral sequence (as exem- 

O I plified e. g. by empirical nearest-neighbor spacings histograms of various quantum systems). The 

^ ■ pertinent Markov diffusion-type processes (with values in the space of spacings) share a general 

form of forward drifts b{x) — -^^^^ — x, where a; > stands for the spacing value. Here A'^ = 2, 3, 5 
^Nj , correspond to the familiar (generic) random-matrix theory inspired cases, based on the exploita- 

^ ' tion of the Wigner surmise (usually regarded as an approximate formula). A'^ = 4 corresponds to 

^>0 ' the (non-generic) non-Hermitian Ginibre ensemble. The result appears to be exact in the context 

\l , of 2 X 2 random matrices and indicates a potential validity of other non-generic N > 5 level 

' repulsion laws. 

m 
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PACS Numbers: 02.50,03.65,05.45 



^ ; 1 Regular versus irregular in quantum chaos 

o 

^ ' The vague notion of so-called quantum chaos, normally arising in conjunction with semiclassical quan- 

tum mechanics of chaotic dynamical systems p| , currently stands for a key- word capturing continued 
5-H ' efforts to give a proper account to what extent quantization destroys, preserves, or qualitatively re- 

produces major features of classical chaos. There is no general agreement about what actually is to be 
interpreted as "quantum chaos" or its definite manifestations. That in part derives from an inherent 
ambiguity of quantization schemes for nonlinear, possibly nonconservative, driven and damped classi- 
cal problems, and is intrinsically entangled with a delicate reverse problem of a reliable (semi) classical 
limit procedure for once given quantum system. Other origins of this elusiveness seem to be rooted 

'Presented at the XIV Marian Smoluchowski Symposium on Statistical Physics, Zakopane, Poland, September 9-14, 
2001 
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in the diversity of meanings attributed in the hterature to the concepts of regular and/or irregular 
behaviour of a physical system, irrespective of whether it is classical or quantum. 

Mathematical definitions of classical chaos emphasize an apparent appearance of randomness in 
deterministic situations, j2| ^. That involves a deep question of when specific features of a physical 
system appear to be (or can be interpreted as) random. Quite typically, while dealing with an irregular 
behaviour, we need to quantify an interplay between chance and order in terms of suitable measures 
of randomness (if random, then "how much"?), @, H, 0]. 

Disorder, irregularity and randomness are casually perceived as synonyms and are interpreted to 
stay at variance with notions of order and regularity. Albeit order and randomness may as well coexist 
as " two faces of the same mysterious coin" , |^ . 

One of basic problems in the quantum chaos theory is to establish whether the classical order - 
disorder interplay induces any unambiguous imprints ("signatures of chaos", [^) in quantum systems. 
In this context, familiar concepts of regular and irregular spectra [||, |^, |l^ were coined to characterize 
distinctive differences between semiclassical distributions of energy eigenvalues for generic quantum 
systems. The term "generic" basically means "more or less typical", since one excludes from consid- 
erations all systems which do not behave properly, although there are many of them. (More stringent 
definition invokes suitable symmetry properties of the quantum system.) 

Nowadays it is clear that an irregularity alone of any particular spectral series (possibly interpreted 
in terms of an irregular sequence of consecutive enery or quasi-energy levels) is not an adequate 
criterion for quantum manifestations of chaos. In fact, semiclassical spectra corresponding to many 
1^ classical systems, be them integrable or chaotic, have an irregular appearance. That was the 
motivation for attempts to classify such spectra in terms of the " degree of randomness of the sequence 
of eigenvalues" , |11 , cf. also [|| and for related argumentation. 

Let us however stress that a primordial question of whether a given energy level sequence can at 
all be regarded as random has been left untouched. Seldom one may have at disposal a complete 
analytic information about quantum spectra. Usually some experimentation is needed to extract the 
data and most of available spectral information relevant for quantum chaos studies, comes either from 
a genuine experiments (microwave analogs of quantum billiards, realistic nuclear data) or computer 
simulations, always with a definite beginning and an end. In each case that produces a finite string 
of data and it is known that no finite sequence can be interpreted as truly random. Fortunately, if a 
data sequence generated by a stochastic process of any origin (deterministic algorithms included) is 
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sufficiently long, then it will always satisfy a test for randomness with fine-tuned confidence level, cf. 

We emphasize an issue of randomness, because various probability laws (and densities of invariant 
measures) are omnipresent in the quantum chaos research. In view of that, a stochastic modeling will 
be our major tool in below. 

A rich class of classically integrable (hence regarded as regular) systems displays random -looking, 
locally uncorrelated sequences of energy eigenvalues, [|[ |l4[ However there are well known classi- 
cally chaotic (hence regarded as irregular) systems whose quantum spectral statistics appears to bear 
no distinctive imprints of classical chaos and look appropriate for the completely integrable case, iQ . 
Therefore, a supplementary rule is necessary to typify various classes of spectral irregularities and of 
the involved types of randomness (being random, but possibly "random otherwise"), if those are to 
be interpreted as consequences of irregular characteristics of the related classical system. 

A possible hint might have originated from discriminating between the spatial regularity and 
irregularity of the corresponding quantum eigenfunctions. It is the spatial pattern of wave functions 
that appears to have a decisive effect on the spectral pattern of eigenvalues, jl^, However, a 
minor obstacle still persists: not all classically ergodic systems (irregular case of Ref. |^) would 
semiclassically yield irregular eigenfunctions, |lO| . Consequently, one usually tries to narrow the 
class of quantum system that are suspected to show undoubtful " signatures of chaos" to those which 
have irregular eigenfunctions, with no specific reference to their classical (chaotic or non-chaotic) 
behaviour. In this class one ultimately attempts to specify those systems which remain in a consistent 
semiclassical relationship with their chaotic classical partners. Those systems quite justifiably would 
deserve to be named generic and would more or less naturally fall into various spectral universality 
classes, in accordance with the random-matrix classification scheme, |l8[ |l9|. 

Under rather plausible assumptions, | |To| , quantum systems with spatially irregular wave functions 
were found to exhibit level repulsion, hence to "avoid" degeneracies which is basically an indication 
of non-integrability, hence not necessarily that of any links with chaos. 

We recall that an opposite spectral effect of level clustering, combined with the conjectured absence 
of correlations between levels, is characteristic for a large class ("almost all" according to 



see however 20 ) of classically integrable systems. Typically they display so-called Poisson statistics 
(strictly speaking, there is an exponential law of probability involved [pT|) of adjacent level spacings: 
small spacings are predominant and there is enough room for multiply degenerate levels. One says 
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then | |22| that energy levels occur in a completely random way via a Poisson process on the energy 
axis. 

Therefore, level repulsion, when regarded as an emergent spectral symptom of level correlations 
(usually interpreted as a certain departure from purely random behaviour), may be viewed as a 
necessary condition to deal either with quantum imprints of classical chaos or, in the least, with a 
classically nonintegrable phase-space irregularity. 

Interestingly enough, this viewpoint finds some support in the discovery of pseudointegrable sys- 
tems (variously shaped billiards, sometimes with singular scattering obstacles) which appear to be 
neither integrable nor chaotic, but give rise to various forms of "wave chaos" while quantized, ||2^, p^ . 
The corresponding distribution of adjacent level spacings is named semi-Poisson and combines vari- 
ous forms (including fractional powers) of level repulsion with Poisson (exponential) statistics, hence 
purely random behaviour in the spectral series. The repulsion phenomenon is here a manifestation of 
the topologically complicated phase space (an invariant manifold is not topologically equivalent to a 
torus but to a higher genus manifold), which was conjectured to preclude integrability and thus the 
standard torus (EBK) semiclassical quantization, see e.g. p^-pGj. 

In view of the wide usage of such terms like "universality" in the quantum chaos literature, one 
should always keep in mind that harmonic oscillators display level repulsion, although they seem 
to be exemplary cases of classical and quantum regularity at their extreme. Another spectacular 
exception is the hydrogen atom spectrum. Like all higher dimensional harmonic oscillators, or a square 
billiard |2^, |l^ , the Coulomb spectrum belongs to a distinctive group of " pathologically nongeneric" 
spectral problems, 

Nonetheless we shall confine our attention to the suggestive, random matrix theory universality 
classification that is considered to be faithful for local fluctuations in quantum spectra of (generic) 
systems that display global chaos in their classical phase spaces. Our hunch is to mimic (or rather 
extract) those features of the level-spacing classification which may bear imprints of pure randomness 
or in reverse - depart from randomness. 

Studying classical manifestations of chaos in terms of probability measures (including their densities 



or distributions and their dynamics) is a respectable strategy, |31 . In quantum theory, in view of 
Born's statistical interpretation postulate alone, probability measures are ubiquitous. On the other 
hand, various probability laws and distributions naturally pervade the familiar random-matrix theory 
[l8|] . This statistical theory of spectra, models a symmetry -limited spectral disorder in terms of 
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statistical ensembles of complex quantum systems (e.g. heavy nuclei). Apart from an ensemble input, 
random-matrix theory forms a convenient vehicle to interpret spectroscopic properties of a concrete 
(single !) quantized version of a complex classical model. (The classical complexity notion refers e.g. 
to the phase-space organization specific to a system and various complication degrees of its dynamics 
related to ergodicity, mixing and exactness.) 

However, one should keep in mind that the universality hypothesis in the context of quantum chaos 
proper, derives from exploiting a spectral afhnity of an ensemble of large (with size ultimately growing 
to infinity) random matrices, with a once given individual quantum system (take into consideration 
the Sinai billiard or (a)periodically kicked pendulum/rotator). Therefore, we may justifiably ask how 
an individual (Hamiltonian or Floquet-type) quantum eigenvalue problem may capture all conceivable 
statistical properties of suitable random-matrix ensemble spectra ? Told otherwise, how may we justify 
a comparison of a statistical ensemble of disordered spectral series with one only specific energy (or 
quasi-energy) level sequence of an a priori chosen quantum system ? 

To our knowledge this conceptual obstacle, except for preliminary investigations of Ref. has 
not received much attention in the quantum chaos literature. A partial answer to that question, 
points towards certain ergodicity properties appropriate for models of the parametric level dynamics 
(Coulomb gas, plasma or else, evolving in "fictitious time" ), that provide a reinterpretation of random- 
matrix theory in terms of an equilibrium statistical mechanics for a fictitious iV-particle system (with 
N allowed to grow indefinitely). 

In the framework of random-matrix theory, an ergodic problem for Gaussian ensembles was ana- 
lyzed long ago in Ref. ||29|| , with a focus on the ergodic behaviour for the eigenvalue density and fc-point 
correlation functions of individual random matrices and their statistical ensembles. That involves a 
local version of the ergodic theorem, where e.g. the spectral averaging over a finite energy span of the 
level density is compared with the matrix ensemble mean of the level density. That suggests analogies 
with disordered many-body quantum systems where ensemble averaging is a standard analytic tool, 
while for an individual system, only an energy averaging should be employed, [ pO[ . 

We know that the distribution of spacings of highly excited quantum systems may involve definite 
laws like e.g. the exponential or Wigner-type distributions. Such laws may be related to definite 
stochastic processes as invariant measures, in particular as asymptotic measures to which the process 
does relax. 

We shall focus upon the parametric dynamics (parametric interpolation) scenario for the nearest- 
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neighbor spacing distributions of irregular quantum systems where asymptotic invariant (with respect 
to the parametric process) probabihty measures are ultimately involved. (Let us mention that a 
concept of parametric dynamics involves the possibly troublesome "fictitious time" parameter. Its 
possible interpretation is that of a running coupling constant measuring the strength of the chaotizing 
perturbation, or more generally that of a "complexity parameter" whose growth to infinity gives 
account of the complexity increase in the spectral properties of a quantum system. 

Disregarding the origins of randomness in diverse settings, we shall take the view that stochastic 
processes are mathematically appropriate models when the time evolution (parametric "dynamics" 
being included) of random phenomena is involved. Whenever probability laws are in usage, random 
phenomena and stochastic processes are always at hand, |3l|, |3^ . 

The major difference of our strategy, if compared to other approaches, amounts to considering 
exclusively the parametric evolution (relaxation) towards equilibiria of nearest neighbor spacing dis- 
tributions as the major source of probabilistic information. We arrive here at prototype invariant 
measures of limiting stationary stochastic processes. We do not invoke any explicit eigenvalue (e.g. 
a solution of the spectral problem for the quantum system or the related random-matrix model) nor 
level dynamics input, since those data prove to be irrelevant for the primary ergodic behaviour that 
is displayed by the adjacent spacing probability densities. At least in the considered approximatio 
regime, where probability densities surmised by Wigner are assumed to be adequate (in reality, they 
have the status of reliable approximate formulas). 

An exploitation of ergodicity (in fact strong mixing and/or exactness, Q) of certain 

(parametrically evolving) Markovian stochastic processes is here found to provide a supplementary 
(probabilistic) characterization of quantum signatures of chaos. 

2 Poissonian level sequences 

2.1 Exponential random variable and semi-Poisson laws 

"Poissonian" matrix ensembles with independent random diagonal elements are often used to model 
spectral properties of integrable Hamiltonian systems (we disregard an issue of various, even quite 
remarkable, deviations from an exact Poisson-type statistics, |2^, 0). Indeed, many regular (inte- 
grable) systems, semiclassically exhibit spacings between adjacent energy levels which are distributed 
according to the exponential probability density p{s) — expi^—s) on where we tacitly assume a 
normalization of the first moment (mean spacing) of the probability measure (hence the unfolding of 
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the energy spectrum, ^, |T^). 

A canonical statement in this respect, conveys a message that "for generic regular systems" 
p{s) is "characteristic of a Poisson process with levels distributed at random" and "the levels are not 
correlated" . (A discussion of serious violations of the Berry - Tabor conjecture can be found in Ref. 

Since the regular spectrum is perfectly deterministic and for each set of quantum numbers the 
corresponding energy level is obtained from an explicit formula (via Einstein-Brillouin-Keller semi- 
classical argument, or directly by solving the spectral problem for e.g. rectangular billiard) it is far 
form obvious that levels may come as a realization of a random variable. Even though probability 
distributions are thought to arise in near classical quantum systems, when the number of levels in any 
range of energy is very large (and indefinitely increases when the classical limit is approached). 

Following Refs. ^ , 16 , let us consider a sequence of numbers (we keep an explicit energy notation. 



although an unfolded sequence is rescaled to be non-dimensional): 

i 

E,+i ^ E,+x^ = Eo + Y^ xj (1) 

where -Eq = and a;j withj = 1, 2, .. are outcomes of independent trials of the exponentially distributed 
random variable X taking values in 

The resulting sequence (i?i, E2, ■■■) of nonnegative numbers is a particular model realization (sam- 
ple) for what is commonly named a Poisson spectrum. Here, randomly sampled (independent, in 
accordance with the exponential distribution) increments Xi — Ei^i — Ei play the role of adjacent 
level spacings. Let us emphasize that the Poissonian random-matrix ensemble would comprise all 
possible sequences of the above form, each obtained as a result of independent sampling procedures. 

At this point let us turn to an explicit probabilistic lore (cf. ^) whose absence is painfully 
conspicuous in major quantum chaos publications. 

Let Xi,X2, ... be independent random variables with common for all exponential probability law 
fj.{x) = aexp{~ax), a > with mean ^ and variance 

Furthermore let us denote Sn = Xi + X2 + ... + A"„, n = 1, 2, .... Then the random variable S'„ 
has a probability density: 

a" 

Pn{x) = —exp{-ax) (2) 

(n- 1)! 

coming from an (n-l)-fold convolution of exponential probability densities on The law (2) is 
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infinitely divisibile, |2^, 

Pn+m{x) (Pn * Pni){x) = I Pn{x - y)pm(y)dy (3) 



where pi{x) — and n, to = 1, 2, .... 

In particular, note that Xi + Xj for any e N has a probability density 

P2{x) — X exp{~ax) . (4) 

which upon setting a = 2 and x = s stands for an example of a semi-Poisson law P{s) = 4s exp(—2s), 
1^, which has been identified to to govern the adjacent level statistics for a subclass of pseudoin- 
tegrable systems. 

It is also obvious that other (plasma-model related, [^) semi-Poisson laws come directly from 
distributions appropriate for Sn- For example, S3 has a density ps^x) which upon substituting a = 3 
and x ~ s gives rise to P{s) = exp{—3s). Analogously, S5 yields p^ix) and upon setting a ~ 5 
implies P{s) = ^s-^expi-bs), cf. Eq. (36) in Ref. |2§. 

2.2 Gaussian regime 

Both in the c^uantum chaos and random-matrix theory contexts, the regime of n >> 1 is of utmost 
importance. Since the primary random variable X has an exponential density with mean A* = ^ ^^^d 
variance cr^ = we stay within the conditions of the central limit theorem, First of all we 

know that for every e > 0: 

P[\-S,,-fi\>e]^0 (5) 
n 

when n ^ 00. Hence -^Sn ^ fJ- with probability 1. 
Furthermore, we have: 



P[a < ^^-^ <b]^ ^= / expi- ^""^ f dx (6) 



To give a pedestrian intuition about the above formal observations, let us ask for a probability 
that there holds 

for any a > 0. In the regime of large n, an integral exp{—y^) dy gives a reliable answer. The 

same integral determines the probability that \Sn — 'n^\ < aay/n, hence tells us how 5„ fluctuates 
about n/i (and -^Sn about /i) with the growth of n. 
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2.3 Whence Poisson process on the energy axis ? 

The probability density of the random variable Sn allows us to evaluate a probability that the n-th 
level energy value En is actually located in an interval [E, E + AE] about a fixed nonnegative number 
E. It is easily obtained by redefining the previous Pn{x), cf. Eq. (2): 

P[E<Sn<E + AE] = _ exp{~aE) AE = (8) 

[n 1 j. 

^ Texp{~S)AS = PniS)AS 



(n-iy. 

where x — E, S = ^-f^ and ^ =< E > is the mean adjacent level spacing. The probability density 
Pn{S), in Ref. is interpreted as "probability density for finding the n-th neighbor of a level in the 
distance increment [S, S + dS] , for a stationary Poisson process" , while in Ref. , while denoted 
E{k, L) — ^ E{n — 1, S) where L is replaced by our S, stands for the "probability that inside an interval 
of length S we find exactly rt — 1 levels" . 

Since E(k, L) has the form of a standard Poisson probability law with mean- value and variance L, one 
may also follow to tell that "if they are on the average L events, then the probability to actually 
observe k events is given by Epoisson{k, L) — j^exp(~L)" . 

Indeed, if E' > is a fixed energy value and we ask for a probability that there are exactly n energy 
levels below E, then probability distributions for Sn and Sn+i combine together to yield the Poisson 
distribution with mean aE: 

P[Sn <E< Sn+i] = P[N{E) =n] = &lexp{-aE) (9) 

n\ 

In this connection, let us recall that a random variable N taking discrete integer values 0, 1, 2, ... 
is said to have Poisson distribution with the mean (and variance) A if the probability of N = k reads 
P[N = /c] ^exp{-\). Clearly, J2T=o ^ k] ^ 1 and E[N] = J^^ kP[N = fc] = A. 

Let us however stress that no explicit Poisson process was involved anywhere in the above, since 
its precise mathematical definition |2^, ^ refers to a counting process with a one parameter family 
of random variables [Nt = N{t) = n] = [Sn < t < Sn+i] obeying the Poisson probability law for all 
t e R+: 

P[Sn < t < Sn+i] - P[Nt - n] - L^exp{-at) . (10) 

ni 

The Poisson process has stationary independent increments: TVt^ , Nt^ — Nt-^ , ... for < ii < t2 < ■■■■ 
with the Poisson probability distribution for each increment: 

P[Nt - Ns = n] = \^^i—^exp[-a{t - s)] (11) 



where Nq = with probabihty 1. Here, by denoting P„(t) = P[Nt = n] and P„(i— s) = P[Nt — Ns ~ n] 
we easily check that /J Pn{t — s)P,n{s){ads) = Pn+,n{t)- 

The related intensity (parameter, mean) of the Poisson "process" equals E[Nt] = at and displays 
the linear growth when t increases. Notice also that — s- a with probability 1 as t — s- oo. 
(The Poisson process is a particular example of a Markovian process in law , |Q . We deal here with a 
temporally homogeneous process iVj, t > associated with an infinitely divisible probability distribu- 
tion fJ,{k) — j^exp{—c), The process in law is here recovered by simply setting /i*(fc) = ^^^j* exp{—ct) 
where fi^{k) — /i(fc).) 

The Poisson process involves time dependent probabilities: Po{t) = exp(—at), Pi{t) — at exp{—at), 
which should be compared with previous outcomes for the exponential random variable. By re- 
calling Eq. (2) we immediately arrive at a formal identification of probability distributions: 

Pk+iit) ^aPk{t). (12) 

In the above, the exponential probability desnity is labelled by time t. Let us stress that pk+i{t) At 
stands for a probability that the random variable Sk+i takes its value in the interval [t, t + At], while 
Pfe(t) is a probability that Nt = n. 

Notwithstanding, Eq. (9) is formally identical with Eq. (10), and therefore we can in principle 
vary the parameter E, so setting (9) in a direct equivalence with a parametric (evolving in fictitious 
time) Poisson process. This formal equivalence underlies a Poisson process lore of the quantum chaos 
literature. 

Instead of paying attention to the exponential probability rule which is responsible for the ran- 
domness of the collection of " time" instants on i?+ , one is tempted to tell that it is the Poisson process 
which dictates those rules of the game. The standard way of thinking refers to the "observation of 
the number of signals recorded up to an instant t (actually, number of jumps of Nt or the number of 
levels that are below E) , . 

2.4 Ergodicity 

Sample paths of the Poisson process Nt are nondecreasing functions of t with integer values. If we 
attempt to draw a sample path, we begin from the value Nt = which is maintained up to the time 
instant Si = ti when the jump occurs to Nt = 1. This value stays constant up to the time 5*2 = xi+X2- 
Then, a new jump to A^f = 3 occurs, and that value survives until S3 ~ xi + X2 + X3 is sampled. The 
sample path construction for the Poisson process strictly parallels a time series construction in terms 
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of points on i?"*" at which jumps of Nt occur. Intervals between consecutive time instants form the 
sequence {xi, X2, x^...) of adjacent level spacings. 

On the other hand, it is Eq. (1) which provides us with a concrete sample sequence of levels 
i?2, ■••); drawn in accordance with the exponential probability law for adjacent level spacings 
Xi,i > 1. Thus , the set of all realizations of the random variable E = (Si, S2, ■■■) comprises a 
statistical ensemble of sample sequences uj : E{u)) = {Ei, E2, ■■■)■ In fact, those sequences exemplify 
the Poissonian ensemble of spectra. 

(If we set a ~ I, then a connection with the standard Poissonian reasoning in the random-matrix 
approach to quantum chaos is immediate. A catalogue of various statistical measures for the Poissonian 
spectra can be found in Jigf.) 

If we would construct a histogram of adjacent level spacings for a single sequence (i?i, E2, , ,) which 
was compiled in accordance with the exponential distribution, the familiar Poissonian shape would be 
revealed. 

As well, the very same picture would emerge if we would randomly collect and make a statistical 
analysis of various finite strings of neighboring energy levels, like in case of the so called nuclear data 
ensemble composition (there e.g. one makes a compilation of 1407 data points from 30 sequences of 
levels experimentally found for 27 different nuclei), ||l9|, p7[ |. 

All that is connected with a primitive at this stage notion of ergodicity of the exponential process . 

Namely, let us consider a one-parameter family {Xn, 71 = 1, 2, ...) of exponential random variables 
as a stochastic process with " discrete time" . Since A"„ are independent random variables with the 
same for all n probability distribution, then for any real function (/ : a; G —>■ ,f{x) G R) such that 
< / >= E[f{Xi)] exists, we have 



for all sample sequences X{lj) ~ (xi,X2, ■■■)■ In that case the random sequence Xn,n > 1 is known 
to be ergodic with respect to f. That is a standard link between the "time average" and "ensemble 
average" , which is here accomplished by means of the exponential probability measure /i. Indeed, as 
often happens in the the context of stationary stochastic processes, ergodicity property allows us to 
replace an average over the set of all realizations of the process at a chosen time instant, by the time 
average evaluated along one sample trajectory. 

If we consider /(X„) ~ Xn for all 71 > 1, then the ergodicity notion refers to limiting properties 
of -^Sn- Accordingly, in view of the law of large numbers (Eq. (5)), Eq. (13) holds true. 




(13) 
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Presently, there is no wonder in the faet that single eigenvalue series of a suitable integrable 
quantum system (like e.g the rectangle billiard of Refs. |l^, |2^, |3^) may be utilized to generate a 
statistical information in (approximate, pO|) afhnity with the ensemble statistics. Numerical research 
involving e.g. about 10® - 10^^ levels for the eigenvalue series = +771^, 7 = 7r/3, cf. [Q, 

allows to generate various statistical data. The nearest neighbor spacing histograms show a very 
close resemblance to the exponential distribution curve, in agreement with the conjecture of Ref. 
Effectively, the eigenvalue sequence of the rectangle billiard can be interpreted as (in fact mimics) a 
sample path E{u}) = (£'i,i?2, ..•) with adjacent spacings Xi distributed according to the exponential 
law. 

A standard (Poissonian) way of thinking in this context, refers to an "observation of the number 
of signals recorded up to an instant t" (actually, jumps of Nt or number of levels that are below E), 
pl| . However, the sample path E{ll>) encodes also a complete information about a sample path of the 
involved exponential process X — (Xi, X2, ■■■)■ 

Our ergodic argument is valid with respect to any chosen sample path X{uj) — (Xi (uj) — xi, X2{lu) = 
X2, ■■■) of X. An ensemble average is provided by xfi{dx) — ^ and that value is to coincide with 
limn^oo-^ irrespective of the particular choice of a sample path uj of the exponential process 

Xn,n> 1. 

Ergodicity property normally embodies the weakest form of complications present in the evolution 
of physical systems, including those modeled by stochastic processes. There is a well established 
catalog of irregular behaviours that the dynamics of any type may exhibit and there are stronger 
types of irregularity than those connected with ergodicity. A corresponding hierarchy of irregularities 
refers to the properties of mixing and exactness Q which will be exploited in below. 

3 Gaussian universality classes: Generalities 

In the random matrix approach we have a priori involved random-looking sequences of energy levels, 
prf , which well agrees with the phenomenology of nuclei where inadequacies of fundamental theoretical 
models are compensated by resorting to statistical matrix ensembles with appropriate symmetries. The 
roots of randomness presumably can be attributed to random deformations of the " shape of the nuclei" 
(bag) in the independent-particle model of nuclear dynamics, | |3^ . An analogue of this reasoning can 
be found in a recent analysis p9| ] of a chaotic system in a cavity (billiard) with a parametric control of 
shape deformations. Then a quantum particle is confined within a continuously deformed boundary , 
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whose parametric dynamics can be as well represented by a stochastic process of any kind. 

A concrete quantum system (like e.g. a spectral problem for concretely shaped billiard) usually 
induces its own unique spectrum and there is no need, nor room for any statistical ensemble of 
systems (unless we shall indeed consider a family of quantum systems with a suitable selection of 
random potentials. We must thus cope with obvious discrepancies underlying otherwise attractive 
affinities (e.g. the universality classes idea for spectral statistics). Useful affinities appear to mask 
quite deep differences between the underlying physical mechanisms. 

It is the level repulsion which is routinely interpreted as a quantum manifestation of classical 
nonintcgrability and ultimately also of chaos, cf. |p^ . 

Normally that is quantified by means of polynomial modifications of the Gaussian probability law 
(in association with the Wigner-Dyson statistics of adjacent level spacings for e.g. unitary, orthog- 
onal and symplectic random matrix ensembles). For completness of the argument, let us list the 



standard formulas: Pi{s) = s^exp{-^), P2(s) = s^^exp{-^) and P^is) = s^^exp{-^), 



corresponding respectively to the GOE, GUE and GSE random-matrix theory predictions. 

Let us point out that for most practical cases the Wigner distributions (albeit exact in the 
2x2 random-matrix case only) are adequate. Typical spacing histograms drawn from experimental 
or numerically generated (quasi)energy spectra are too rugged to allow subtle distictions against the 
n ^ oo random-matrix size related predictions. 

We shall consider mostly the Wigner-type cases, even though neither of those probability laws 
deserves the status of being an exact representation of the real state of affairs. Remember that also in 
the context of random matrix theories the Wigner spacing formulas are approximations that usually 
improve in the large matrix size regime. 

The nearest neighbor spacing distributions, in the random-matrix approach are the secondary 
notions and can be derived from an explicit formula for the joint probability density to find the 
(dimensionless) energy eigenvalues in respective infinitesimal intervals [xi, Xi-\-Axi\ with i — 1,2, N: 



where /3 = 1,2,4 and Cn0 is a normalization constant, g ^7j. The level repulsion has been built 
into the framework from the very beginning and appropriate level spacing distributions (including the 



There were many attempts to provide convincing (and independent from the definite symmetry 
and Gaussian randomness inputs, proper to random- matrix theory)) arguments that would generate 




(14) 



i>j=i 



i=l 



adjacent level case) can be directly evaluated on that basis, ^ pl| . 
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level repulsion through well defined dynamical mechanisms (like e.g. the parametric level dynamics) 
and would lead to statistical predictions as well. A suitable level dynamics scenario may as well 
give rise to the so-called intermediate statistics and possibly a continuous (parametric) interpolation 
among them. 

In the random matrix theory context a radical probabilistic attempt due to Dyson explicitly 
involves the (parametric) Brownian motion assumption for each energy level separately, 

More satisfactory results were obtained by resorting to a fictituous gas of interacting particle 
representatives of individual energy levels. A corresponding many-particle system is then investigated 
at suitable "thermal equilibrium" conditions. Then, without introducing a priori statistical ensembles 
of random matrices, level distribution functions are derived by means of ordinary statistical mechanics 
methods. That approach explictly involves the many-body Hamiltonian (Calogero model): 

whose squared ground state function (equlibrium measure density) has the form (14), 

Apart from that, explicit quantum mechanical investigations for billiard- type systems provide hints 

about the potential importance of interpolation studies, especially since various intermediate types of 

statistics were reported to occur, see e.g. [|l[ 

There are two basic approaches to an interpolation issue. One refers explicitly to random matrix 

theories and their "affinity" with quantimr chaotic systems, Another refers to the fictitious 

gas, interacting many-body analogy, ||l^, ^ 48, |^. Recently, a related short-range plasma 



model was proposed to analyze an emergence of the "pseudo-Poisson statistics", p6{ . 

4 Parametric dynamics of adjacent level spacings 
4.1 Markov processes defined through their invariant measures 

Once we have encountered probability densities on the positive half- line in R^, it is rather natural to 
investigate a general issue of parametric stochastic processes which would provide a dynamical model 
of level repulsion in an irregular quantum system and generate at the same time spacing densities 
as those of asymptotic invariant (equilibrium) probability measures. Such random processes clearly 
must run with respect to the previously mentioned " fictitious" time-parameter and take values in the 
set of all level spacings which are appropriate for a complex quantum system or the corresponding 
random-matrix ensemble. 
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Effectively, we wisli to introduce a Marlcovian diffusion-type process wliicli miglit stand for a 
reliable approximation of a random walk over level spacing sizes. 

For future reference let us mention that in the regime of equilibrium (when an invariant measure 
appears in the large " time" asymptotic) , a sample path of such random walk would take the form of an 
ordered sequence of spacings which are sampled (drawn) according to the prescribed invariant prob- 
ability distribution. That is precisely one explicit example of the ladder of energy levels, understood 
as a random sample drawn from a suitable ensemble. 

An analysis of statistical features of this spectral sequence involves an ergodicity notion to stay in 
conformity with the ensemble evaluation of various averages (carried out with respect to the invariant 
density), |§. 

We shall consider the previously listed GOE, GUE and GSE probability densities on i?+ (up 
to suitable rescalings !) as, distorted in view of the spacing size normalization, asymptotic invari- 
ant densities of certain parametric Markovian stochastic processes whose uniqueness status can be 
unambiguously settled. 

Let is begin from the observation that probability densities on of the characteristic form f{x) ^ 
xexp{—^), g{x) ~ x'^exp{—^) and h{x) ~ ^exp{—x'^) appear notoriously in various quantum 
mechanical contexts (harmonic oscillator or centrifugal- harmonic eigenvalue problems), cf. p3| , ^ 
K3, M, M. Notwithstanding, as notoriously they can be identified in connection with special classes 



of stationary Markovian diffusion processes on |57 



Anticipating further discussion, let us consider a Fokker-Planck equation on the positive half-line 
in the form: 

dtp^^Ap~V[^-x)p] (16) 

which may be set in correspondence with the stochastic differential equation dXt ~ {^-^ — Xt)dt+dWt 
formally valid for a random variable Xt with values contained in (0, oo). Here /3 > and Wt represents 
the Wiener process. 

Accordingly, if pq{x) with x E is regarded as the density of distribution of Xq then for each 
t > the function p{x, t), solving Eq. (16), is the density of Xf. In view of a singularity of the forward 
drift at the origin, we refrain from looking for strong solutions of the above stochastic differential 
equation and confine attention to weak solutions only and the associated tractable parabolic problem 
(16q) with suitable boundary data, cf. pTf . 

In all those cases a mechanism of repulsion is modeled by the - term in the forward drift expression. 
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The compensating harmonic attraction which is modeled by the —x term, saturates the long distance 
effects of repulsion-induced scattering and ultimately yields asymptotic steady (stationary) probability 
densities. 

To interpret a density p{x) as an asymptotic (invariant) density of a well defined Markovian 
diffusion process we shall utilize the rudiments of so-called Schrodingcr boundary and stochastic 
interpolation problem, ^ see also when specialised to invariant measures. 

Let us notice that both in case of the standard Ornstein-Uhlenbeck process and its Bessel (radial) 
variant, we have emphasized the role of a stochastic process with an asymptotic invariant density. To 
deduce such processes, in principle we can start from an invariant density and address an easier issue 
of the associated measure preserving stochastic dynamics. Next we can consider whether the obtained 
process would drive a given initial density towards a prescribed invariant measure (in that case we 
can tell about an asymptotic state of equilibrium to which the process relaxes) . That feature involves 
the notion of exactness of the related stochastic process, whose straightforward consequence are the 
properties of mixing and ergodicity of the corresponding random dynamics, pi] ]. 

There is a general formula ||5^, |l], ^ relating the forward drift of the sought for stationary process 
with an explicit functional form of an invariant probability density. We confine our attention to Markov 
diffusion processes with a constant diffusion coefficient, denoted D > 0. Then, the pertinent formula 
reads: 

b{x) = 2D-^. (17) 

In particular, for the familiar Ornstein-Uhlenbeck process we have p^^^{x) — (i)^/^ea;p(— ^) and 
Z) = ^, so we clearly arrive at b{x) = ~x as should be. Quite analogously, in case of the GUE-type 
spacing density, we have D = ^ and p^/'^{x) — ^j^x exp(— ^). Thus, accordingly b{x) — ^ ~ x. 

The very same strategy allows us to identify a forward drift of the Markovian diffusion process 
supported by the GOE-type spacing density. By employing p^/^{x) — \/2xexp{—^) and setting 



D — ^ we arrive at the formula: b{x, t) — 



X. 



We immediately identify the above forward drifts with the ones appropriate for the time homo- 
geneous radial Ornstein-Uhlenbeck processes, with a corresponding family of {N > 1 and otherwise 
arbitrary integer) transition probability densities, p7[ : 

Pt{y, x) = p{y, 0, X, t) = 2x^-1 expi-x^y (18) 

1 {x^ + y^)exp{-2t) 

1 - exp(-2i) l-cxp(-2i) ^' 
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r / ,M-a. , 2xyexp{-t) 

'"^ l-exp(-2t) ^ 

where a — ^^"^ and Ia{z) is a modified Bessel function of order a: 

^-(-')-g (M)r(V+o^i) 

while the Euler gamma function has a standard form T{x) = cxp{~t)t^~^dt. We remember that 
r{n + l)=n\ and r(l/2) = V^. 

The resultant forward drift has the general form: 

iV - 1 

b{x)^^^-x (20) 

and corresponds to f3 — N ~ I. 

By setting N — 2, and then employing the series representation of Io{z), we easily recover the 
asymptotic invariant density for the process: limt-^ocP{y,0,x,t) — 2xexp(— x^). 

We can also analyze the large time asymptotic of p{y, 0, x, t), in case of = 3 which gives rise to 
an invariant density in the form: ■^x'^exp[—x^). That obviously corresponds to the GUE-type case 
with b[x) = ^ — X. 

When passing to the GSE case, we are interested in the Markovian diffusion process which is 
supported by an invariant probability density p{x) — ^{3/2) cxp(— a:^). Let us evaluate the forward 
drift of the sought for process (we set D — b{x,t) = | — a;. Clearly, we deal here with a 
radial Ornstein-Uhlenbeck process corresponding to — 5. The transition probability density of 
the process displays an expected asymptotic: /imf^oop(y, 0, x, i) = -^x* exp(— x^). Here we have 
exploited r(l/2) = 0r to evaluate r(3/2) = i0r. 

The above formulas allow us to formulate a hypothesis that further non-generic repulsion laws may 
be appropriate for quantifying quantum chaos. Straightforwardly, one can verify that our transition 
probability densities refer to asymptotic invariant densities of the form: 

p{x) = ^:^^x''-'exp{-x'). (21) 

In particular we get a direct evidence in favor of iV = 4, i. e. b{x) = ^ ^ 2;, universality class 
which in fact corresponds to the Ginibre ensemble of of non-Hermitian random matrices, where a 
cubic level repulsion appears: p{x) = 2a;'^ exp(— x^) (this formula is exact for 2x2 random matrices). 

In principle, processes corresponding to any N > 5 may be realizable as well, and thus the related 
higher-power level repulsion might have relevance in the realm of quantum chaos. 
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In all considered cases, an asymptotic invariance of probability measures (densities) is sufficient to 
yield ergodic behaviour. For each value of iV > 1 we deal with an independent repulsion mechanism, 
albeit all of them belong to the radial Ornstein-Uhlcnbeck family. 

We have thus identified a universal stochastic law (in fact, a family of the like) behind the functional 
form of basic, Wigner surmise inspired, spacing probability densities appropriate for quantum chaos. 

Let us emphasize at this point that one should keep in mind a number of possible reservations 
coming from the fact that neither of "universal" or "generic" laws can be regarded as a faithful 
representation of a real state of affairs. Usually exact laws are derived for two by two (hence of the 
small size !) random matrices, and are known to reappear again (at least in the generic cases) as 
approximate spacing formulas in the large random-matrix size regime. That in turn allows to achieve 
a correspondence with semiclassical quantum spectra of complex systems. 

There is no obvious explanation of a physical meaning of the integer parameter N in the radial 
stochastic process scenario. One hypothesis comes from the random-matrix theory, where f3 = N —I = 
1,2,4 would correspond to a number of independent components of a typical matrix entry which is 
decided by the underlying symmetry of the problem (GOE, GUE, GSE). That can be presumably be 
extended to the case of iV = 4 and possibly all < 8. 

4.2 Link with Calogero Hamiltonian 

Previously we have indicated that a common mathematical basis for various level repulsion mechanisms 
appropriate to quantum chaos is set by the Calogero-Moser Hamiltonian. At the first glance, our 
stochastic arguments may leave an impression that something completely divorced from that setting 
has been obtained in the present paper. However things look otherwise and our theoretical framework 
proves to be compatible with standard techniques for spectral analysis of complex quantum systems. 

It is peculiar to the general arguments of Refs. ^ that invariant probability densities give 
rise to measure preserving stochastic processes in a fully controlled way. One of basic ingredients 
of the formalism is a proper choice of Feynman-Kac kernel functions, which are the building block 
for the construction of transition probability densities of the pertinent Markov processes. Feynman- 
Kac semigroup operators (and their kernels) explicitly involve one particle Hamiltonian operators as 
generators (in less technical terms one may think at this point about rather standard transformation 
from the Fokkcr-Planck operator to the associated self-adjoint one, |^). 

For stationary processes, a general formula relating forward drifts b{x) of the stochastic process 
with potentials of the conservative Hamiltonian system reads (we choose a diffusion coefficient to be 
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equal i), g, [58 



V{x)^^{b'+V-b). (22) 



Upon substituting the general expression for b{x) we arrive at: 

V{^)-h^^-^~{P + l) + -'] (23) 



where [3 = N — 1. This potential function enters a standard definition of the one particle Hamiltonian 
operator (physical parameters have been scaled away): 

H=-\li + V{x) (24) 

where A = -j^. The operator (24) with V{x) defined by (23) is an equivalent form of a two-particle 
(actually two-level) version of the Calogero-Moser Hamiltonian, cf. fs^ . 
Indeed, the classic Calogero-type problem is defined by 

with the well known spectral solution. The eigenvalues read E„{f3) — 2n + 1 + i[l + /3(/3 — 2)]^/^, 
where n > and f3 > —I. 

By inspection we can check that all previously considered iV = 2, 3, 4, 5 radial processes correspond 
to the Calogero operator of the form H—Eq where Eq is the ground state (n=0) eigenvalue. Its explicit 
foom relies on the choice of /3 and by substituting f3 — 1, 2, 3, 4 we easily check that 

EoW = 1 + i[l + /3(/3 - 2)]i/2 = i(/3 + 1) . (26) 

Accordingly, all considered radial processes arise as the so-called ground state processes associated 
with the Calogero Hamiltonians (squared modulus of the ground state wave function stands for the 
pertinent probability density). Let us recall that the classic Ornstein - Uhlenbeck process can be 
regarded as the ground state process of the harmonic oscillator Hamiltonian operator. That by the 
way corresponds to choosing iV = 1 i.e. /3 = in the above, plus allowing the whole of to the 
process. Like in the standard OU process case, radial OU processes share the property of exactness 
(while driving any initial density towards suitable equilibrium) and hence ergodicity. 

5 Discussion 

Our motivations were essentially probabilistic and spectral series with spacing densities governed by 
Wigner-type laws have emerged in the course of a parametric stochastic process that relaxes towards 
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equilibrium (invariant measure). Such series have thus a definitely random origin. It is clear that an 
approximate value of Wigner densities indicates nonrandom input in realistic cases. 

Let us point out that in the standard matrix-theory framework Dyson' s " threefold way" is based on 
the demonstration that on general (invariance under symmetry) grounds only three basic ensembles 
(orthogonal, unitary and symplectic) matter. Hence, the non-generic repulsion behaviour we have 
discussed before, goes beyond the standard framework (under assumptions of the Dyson theorem the 
non-generic laws are not admissible). Many different ensembles have been used in the literature, but 
their properties were more specific (less general) than the standard GOE, GUE and GSE cases show 
up. 

The spacing distributions we have addressed (Wigner surmise), fail to be correct in general. The 
true random-matrix universal distributions differ from them, albeit the discrepancy is small when 
matrix size is going to infinity, 

In the discussed parametric relaxation process scenario, one may easily implement a transition of 
any initial density towards a concrete asymptotic one with the wealth of intermediate examples (e.g. 
from Poisson to GOE interpolation). In that case, both the initial and terminal distributions refer 
to random seqences of numbers (possible energy eigenvalues) . Ergodicity of the ultimate stationary 
process implies that its sample paths arise as random sequences drawn from the Wigner-type distri- 
bution. Clearly, we cannot expect that such purely random sample sequences would reveal long range 
correlations typical of random-matrix models. 

Acknovifledgement: I would like to thank Karol Zyczkowski for correspondence on semi- Poisson 
laws and for pointing Ref. [po[ to my attention . 
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